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ABSTRACT:- This paper presents a method to convert the Meteorological tabulated data into an 
analytical expression in which it can be used in engineering application. 

The paper uses a Fourier series technique to convert the solar radiation variation during the day for clear 
days and for a beclouded intermittent days to predict an analytical expression for solar radiation at any time 
of that day by finding the weights for both (a n & b n ) parameters for Fourier series upto (n =10). 

It is found an encouraged result that will help the researchers to use Fourier series to convert any tabulated 
data into an analytical expression which gives a useful indication for using such method in solar 
applications. 

The results which are drawn in the plots for (n upto 10) showing the accuracy in coincidence of this series 
increases for the given function as the (n terms) increases from (n = 0 to n =10), for more precision the 
weights for both (a„ & b n ) parameters (n more than 10) should be found. 
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I. INTRODUCTION 

The sun is probably the most important source of renewable energy available today. Solar energy is the 
energy force that sustains life on Earth for all plants, animals, and people. It provides a compelling solution for 
all societies to meet their needs for clean, abundant sources of energy in the future [1]. 

Solar radiation as it passes through the earth's atmosphere undergoes to many factors such as 
absorption, scattering and reflection by air molecules, water vapor and clouds. The direct solar beam that is 
arriving directly to the earth’s surface is considered as direct solar radiation. The total amount of solar radiation 
received and falling on a horizontal surface (which is the summation of the direct solar beam plus diffuse solar 
radiation on a horizontal surface) is called global solar radiation. 

Solar radiation, incident to the earth's surface, is a result of complex interactions of energy between the 
atmosphere and surface [2], 

Accurate measurement devices for solar energy are not available. Devices like Pyranometers are 
capable of measuring solar intensity but, they fail to draw the exact picture of the incident energy. There are 
hundreds of meteorological stations distributed around the world directly or indirectly measuring solar radiation. 
These stations work to derive spatial databases from these measurements with different interpolation techniques 
are used to get the required accurate information or estimation. Other techniques of generating spatial databases 
are solar radiation models integrated within geographical information systems (GIS) [3]. 

In addition, solar energy systems require collecting these records of solar radiation and other 
meteorological data like wind speed, temperature, humidity ratio and cloud which covers the sky. The 
constructed hourly data base where tabulated and arranged in Handbook which are available in computer 
compatible form [4]. 

These tabulated data base is needed from time to time to applied in engineering works, so it is 
necessary to construct an analytical expression to be defined instead of a table of value. Various methods have 
been derived to do this work; the simplest one which is the most obvious in evaluating of the integrals in the 
Euler formulas is done by using the (Trapezoidal Rule). But since for the occasional application that is 
confronting the most researchers and those who are dealing and handling with numerical functions regularly, the 
Fourier expansion is quite satisfactory for such application and it will be discussed in this paper. 

Because so many problems in harmonic analysis involve functions, such as meteorological or 
economic quantities, whose period is either a day or a year, it is customary to assume that data are available at 
interval of (1/12, 1/24, or sometime 1/48 of a period [5]. 
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II. THE RESEARCH GOAL 

The goal of this research is to carry out a study to change the tabulated data base for solar radiation 
intensity into a developed analytical expression or formula that can be used in any engineering applications 
required in solar energy systems or in any meteorological data base such as temperature, wind speed and any 
parameters that take place regularly and in harmonic way. 



III. THEORETICAL ANALYSIS 



Fourier series is that method used to approximate functions with sums of sine and cosine functions. It is 
capable and well suited to analyze periodic and discontinuous functions, such as communication signals and 
alternative current, for solving heat transfer applications, and for many other problems in engineering [6]. 
a 0 n t 2nt nn t nt 2nt 
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We shall consider a function [ f (t) ] of a period (2p) for which values are available at intervals of (1/24): 

2 p p 

At = — = — 

24 12 



Referring to [ figure ( 1 ) ] any function [/(t) ] can be expressed as the sum of an even functions and an odd 
functions, simply by writing: 
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Since [ g (t) ] is clearly even and [ h (t) ] is clearly odd, hence the cosine terms in the expansion of [ f (t) ] are 
just the terms in the half range cosine expression of [ g (t) ], and the sine terms in the expansion of [ f (t) ] in the 
half range sine expansion of [ h (t) ]. 

For the even function [ g (t) ], we have as usual: 
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By applying the expression of trapezoidal rule, with [ At = p /12 ], it is easily can obtain that: 
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But we have: 



Thus: 
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a n = 0.08333/(0) + 0.08333 [/(l) + / (-l)].cos — +.. + 0.08333 [fill) + / (-11)]. cos — 

+ 0.04167[/(12)+ /(-12)] .cosntr ...(7) 

The cosine factors in the last expression can be evaluated once and for all and combined with the other 
numerical factors, including the [ 1/2 ] in the definition of [ g (t) ] to yield a set of weights by which the 
successive values of the sum [/(f) + f (-t) ] are to be multiplied before they are added. In this way, weights of 
[a 0 (n = 0) & aj (n = 1)] for example take the form of: 
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0.08333 /(0) + 0.08333 [f (1) + / (- 1)] + + 0.08333 [/ (11) + / (- 11)] + 0.04167 [f (1) 
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0.08333 /(0) + 0.08049 [/ (1) + / (- 1)] + - 0.08049 [/ (11) + / (- 11)] - 0.04167 [f (1) 
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In similar manner, the other of (a's) weights can be calculated. Table (1) shows the weights involved in 



computation of the first ten of [ a's ]. 



Table ( 1 ): The weights for the first ten of ( a n ). [5] 



Terms 








Weight for determination of ( a 


n) 








ft + ft 


n = 0 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


9o = fo 


0.0833 


0.0833 


0.0833 


0.08333 


0.0833 


0.0833 


0.0833 


0.0833 


0.0833 


0.0833 


0.0833 


3 


3 


3 




3 


3 


3 


3 


3 


3 


3 


fi + fi 


0.0833 


0.0804 


0.0721 


0.05893 


0.0416 


0.0215 


0.0000 


- 


- 


- 


- 


3 


9 


7 




7 


7 


0 


0.0215 


0.0416 


0.0589 


0.0721 


















7 


7 


3 


7 


Z + f-2 


0.0833 


0.0721 


0.0416 


0.00000 


- 


- 


- 


- 


- 


0.0000 


0.0416 




3 


7 


7 




0.0416 


0.0721 


0.0833 


0.0721 


0.0416 


0 


7 












7 


7 


3 


7 


7 






Z + f-3 


0.0833 


0.0589 


0.0000 


- 


- 


- 


0.0000 


0.0589 


0.0833 


0.0589 


0.0000 


3 


3 


0 


0.05893 


0.0833 


0.0589 


0 


3 


3 


3 


0 












3 


3 












f* + f-4 


0.0833 


0.0416 


- 


- 


- 


0.0416 


0.0833 


0.0416 


- 


- 


- 




3 


7 


0.0416 


0.08333 


0.0416 


7 


3 


7 


0.0416 


0.0833 


0.0416 








7 




7 








7 


3 


7 


fs + fs 


0.0833 


0.0215 


- 


- 


0.0416 


0.0804 


0.0000 


- 


- 


0.0589 


0.0721 




3 


7 


0.0721 


0.05893 


7 


9 


0 


0.0804 


0.0416 


3 


7 








7 










9 


7 






fs + f-6 


0.0833 


0.0000 


- 


0.00000 


0.0833 


0.0000 


- 


0.0000 


0.0833 


0.0000 


- 
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3 


0 


0.0833 




3 


0 


0.0833 


0 


3 


0 


0.0833 








3 








3 








3 


fy + f-7 


0.0833 


- 


- 


0.05893 


0.0416 


- 


0.0000 


0.0804 


- 


- 


0.0721 




3 


0.0215 


0.0721 




7 


0.0804 


0 


9 


0.0416 


0.0589 


7 






7 


7 






9 






7 


3 




fa + f-8 


0.0833 


- 


- 


0.08333 


- 


- 


0.0833 


- 


- 


0.0833 


- 




3 


0.0416 


0.0416 




0.0416 


0.0416 


3 


0.0416 


0.0416 


3 


0.0416 






7 


7 




7 


7 




7 


7 




7 


fs + f-9 


0.0833 


- 


0.0000 


0.05893 


- 


0.0589 


0.0000 


- 


0.0833 


- 


0.0000 




3 


0.0589 


0 




0.0833 


3 


0 


0.0589 


3 


0.0589 


0 






3 






3 






3 




3 




fio + f -10 


0.0833 


- 


0.0416 


0.00000 


- 


0.0721 


- 


0.0721 


- 


0.0000 


0.0416 




3 


0.0721 


7 




0.0416 


7 


0.0833 


7 


0.0416 


0 


7 






7 






7 




3 




7 






fll + f-11 


0.0833 


- 


0.0721 


- 


0.0416 


- 


0.0000 


0.0215 


- 


0.0589 


- 




3 


0.0804 


7 


0.05893 


7 


0.0215 


0 


7 


0.0416 


3 


0.0721 






9 








7 






7 




7 


fl2 + f-12 


0.0416 


- 


0.0416 


- 


0.0416 


- 


0.0416 


- 


0.0416 


- 


0.0416 




7 


0.0416 


7 


0.04167 


7 


0.0416 


7 


0.0416 


7 


0.0416 


7 






7 








7 




7 




7 





Similarly, for [ h (t) ], it can easily obtain that: 
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12 



In similar way for [ (n = 0), b 0 will equal zero ] since [ sin ( 0 ) = 0 ], and the weight of ( bi ) for example takes 
the form of: 



hi = 0.08333 [/(l) — /"( — 1)]. sin — — h 0.08333 [/(2) - /(— 2)]. sin llJL + ... ... 

12 6 
1 1 ^ 1 ^ TC 

+ 0.08333 [/ (11)- /(— 11)] . sin — 

b 1 = 0.02157 [/Cl) - /(— 1)] + 0.04167 [/( 2) - /(— 2)] + + 0.02157 [f (11) - /(-ll)] 

The weights involved in computation of the first ten of [ b's ] required for the evaluation of this expression are 
given in table (2). 
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Table ( 2 ): The weights for the first ten of ( b n ^[5] 



Terms 


Weight for determination of [ b n ] 


ft + ft 


n = 1 


n = 2 


n = 3 


n = 4 


n =5 


n = 6 


n = 7 


n = 8 


n =9 


n = 10 


fi - fi 


0.0215 


0.0416 


0.0589 


0.0721 


0.0804 


0.0833 


0.0804 


0.0721 


0.0589 


0.0416 




7 


7 


3 


7 


9 


3 


9 


7 


3 


7 


f2 - f-2 


0.0416 


0.0721 


0.0833 


0.0721 


0.0416 


0.0000 


- 


- 


- 


- 




7 


7 


3 


7 


7 


0 


0.0416 


0.0721 


0.0833 


0.0721 
















7 


7 


3 


7 


f3 - f-3 


0.0589 


0.0833 


0.0589 


0.0000 


- 


- 


- 


0.0000 


0.0589 


0.0833 




3 


3 


3 


0 


0.0589 


0.0833 


0.0589 


0 


3 


3 












3 


3 


3 








f* - u 


0.0721 


0.0721 


0.0000 


- 


- 


0.0000 


0.0721 


0.0721 


0.0000 


- 




7 


7 


0 


0.0721 


0.0721 


0 


7 


7 


0 


0.0721 










7 


7 










7 


fs - f-s 


0.0804 


0.0416 


- 


- 


0.0215 


0.0833 


0.0216 


- 


- 


0.0416 




9 


7 


0.0589 


0.0721 


7 


3 


7 


0.0721 


0.0589 


7 








3 


7 








7 


3 




fe - f-6 


0.0833 


0.0000 


- 


0.0000 


0.0833 


0.0000 


- 


0.0000 


0.0833 


0.0000 




3 


0 


0.0833 


0 


3 


0 


0.0833 


0 


3 


0 








3 








3 








f7 - f-7 


0.0804 


- 


- 


0.0721 


0.0215 


- 


0.0216 


0.0721 


- 


- 




9 


0.0416 


0.0589 


7 


7 


0.0833 


7 


7 


0.0589 


0.0416 






7 


3 






3 






3 


7 


f8 - f-8 


0.0721 


- 


0.0000 


0.0721 


- 


0.0000 


0.0721 


- 


0.0000 


0.0721 




7 


0.0721 


0 


7 


0.0721 


0 


7 


0.0721 


0 


7 






7 






7 






7 






f> - f-9 


0.0589 


- 


0.0589 


0.0000 


- 


0.0833 


- 


0.0000 


0.0589 


- 




3 


0.0833 


3 


0 


0.0589 


3 


0.0589 


0 


3 


0.0833 






3 






3 




3 






3 


fio - fio 


0.0416 


- 


0.0833 


- 


0.0416 


0.0000 


- 


0.0721 


- 


0.0721 




7 


0.0721 


3 


0.0721 


7 


0 


0.0416 


7 


0.0833 


7 






7 




7 






7 




3 




fu - fit 


0.0215 


- 


0.0589 


- 


0.0894 


- 


0.0804 


- 


0.0589 


- 




7 


0.0416 


3 


0.0721 


9 


0.0833 


9 


0.0721 


3 


0.0416 






7 




7 




3 




7 




7 



IV. THEORETICAL APPLICATION 

In this paragraph, and in order to give an example to how apply the Fourier series to serve this method 
a sample of theoretical calculation will put into practice to find the way to convert the tabulated data for solar 
radiation intensity or any other meteorological data into an analytical expression, to do that we should go back 
for any tabulated data for solar radiation intensity. 

We will choose any arbitrary data for clear sky day and the other for partially cloudy day. From the 
Meteorology Database for Baghdad City (2000), the solar radiation intensity of clear sky day for (21 st of 
December) can be tabulated in the following table (3): 



Table ( 3 ): The Solar Radiation for 21 st of December [7], 



Time [hour] 


0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 






















0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


0 


1 


2 


3 


4 




lolar Radiation Intensity [ W/m 2 ] 




0 


0 


0 


0 


0 


0 


0 


3 


1 


2 


3 


3 


3 


3 


2 


1 


3 


0 


0 


0 


0 


0 


0 


0 


0 
















4 


5 


5 


0 


7 


9 


5 


7 


5 


4 


































3 


5 


6 


4 


1 


7 


2 


3 





















The tabulated solar intensity was drawn as shown in [figure (2)], in which [ f (t) & f (- t) ] are 
explained, so the functions [f (t) & f (- 1)] can be calculated easily and tabulated in table (4): 
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» t 



Figure ( 2 ): the solar radiation intensity for (21 st of Dec.) showing [/ ( t ) and / ( - 1 ) ]. 



Table ( 4 ): The values of [ f (t) & f(-t) ]. 



f (-12) 


f(-H) 


f(-10) 


f(“9) 


f(“8) 


f(-7) 


f(~6) 


f(-5) 


f(-4) 


f(-3) 


f(-2) 


f(-l) 


391 


357 


272 


153 


34 


0 


0 


0 


0 


0 


0 


0 



f(0) 


f(l) 


f(2) 


f(3) 


f(4) 


f(5) 


f(6) 


f(7) 


f(8) 


f(9) 


f(10) 


f(U) 


f(12) 


0 


0 


0 


0 


0 


0 


0 


34 


153 


255 


306 


374 


391 



The values of ( a n and b„ ) now can be found by utilizing table (4) in conjugation with equations ( 7 & 9 ), the ( 
a n and b„ ): 

72 72 1172 72 

a n = 0.08333 * 0 + 0.08333 [0 + 0], cos^- + + 0.08333 [374 + 357] . cos — — — 

+ 0.04167 [391+ 391].cosn7r 

b n = 0.08333 [/(l) - /(— l)].sin^-+ 0.08333 [/(2) - /(-2)].sin^- 

12 6 
lln n 

+ 0.08333 [/( 11) - /(-ll)] .sin — 

For ( n = 0 ), [ a 0 ] was found to be equal to ( 194.08 ) and for ( n = 1 ), [ ai ] was found to be equal to ( - 165.71 
)■ 

For ( n = 0 ), [ b„ ] was found to be equal to ( 0 ) since [sin (0) = 0], and for ( n = 1 ), [ bj ] was found to be 
equal to (19.12). 

For other values of ( a„ & b„ ) involved for the calculation for the first ten of them (upto n = 10), they were 
tabulated and shown in table ( 5 ): 



Table ( 5 ): The values of first ten of [a n & b„]. 



Values of ( a n ) 



n = 0 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n= 10 


194.08 


-165.7 


99.18 


-34.03 


-1.41 


6.88 


0.00 


-4.20 


2.83 


0.03 


-1.42 


Values of (b n ) 


n = 0 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


n = 10 


0 


19.12 


-21.67 


7.84 


7.36 


-10.93 


4.25 


3.27 


-4.91 


2.18 


0.42 



Substituting the values the first six of [a n and b„] into equation (1), it yields to an analytical expression for the 
solar radiation [/(t) ] at any required time: 

n t 2n t 6n t n t 

f (t) = 97.04 — a 1 165.7 cos 1 99.18 cos 1 — 0 * cos h 19.12 sin- 

1 1 12 12 

2 n t 6 n t 

— 21.67 sin— — + + 4.25 sin- 



12 



12 



12 



12 
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To obtain for example the value for solar radiation [/(ll)] at time ( t = 11 hour), it can be done simply by 
substituting for ( t = 1 1) in this equation to get its value, which is: 

IT 1 1 "~) TT *11 ^ TT "1 "1 7 T jfc "1 "1 

/ (11) = 97.04- OilbSJ cos + 99.18 cos + + 6.88 cos + 19.12 sin 

7 1 12 12 12 12 

2 n * 11 6 n * 11 

— 21.67 sin 1- + 4.25 sin 

12 12 
, . W 

f (11) = 346.5— r 

m z ^ 

Comparing the actual tabulated value for the solar radiation at ( t = 11 hour ) which is (374 W/nr) with the 
value obtained from the analytical expression 346.5 W/m 2 ), the values are too close for each other. 

V. RESULTS AND DISCUSSION 

In this paragraph, the analysis for this research shows the possibility to predict the solar radiation 
intensity in analytical expression rather than in tabulated data. The analysis can be carried out by choosing the 
solar radiation for a dedicated day which want to convert into analytical expression from a Meteorological 
database and performing the procedure of Fourier series technique to obtain the analytical expression for that 
day. 

For a clear shining day of ( 21 st of December ), the Meteorological database was shown in table (no. 1). 
Performing the Fourier series procedures that was explained in theoretical application paragraph upto ( n = 6 ), it 
shows a good precision coincidence between the tabulated database and the analytical expression. 

Figure (3) shows the solar radiation intensity for (21 st of December) in its analytical expression and 
Meteorological database for (n = 4), the figure shows the coincidence between them for time period of ( t = 6 to 
t = 15 ) is too close, but it is not for the rest of time periods. 

Figure (4) shows the solar radiation intensity for the same day but for other values of (n) for (n = 2 and 
n = 6), it is clear from the figure that as ( n ) increases the coincidence precision increases to obtain a good 
fitness between the values for both of them. 




Figure ( 3 ): the Meteorological data for solar radiation intensity for (21 st of Dec.) and the analytical 

expression for ( n = 4 ). 




Figure ( 4 ): the Meteorological data for solar radiation intensity for (21 st of Dec.) and the analytical 

expression for ( n = 2 & n =6 ). 
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For other day which is not clear or shiny day, and it could be cloudy and sun shining beclouds 
intermittently, it can also using Fourier series to get and analytical expression for that day. 

For instance, the Meteorological data for the solar radiation intensity for (27 lh of Feb.) which is tabulated in 
table (6) can be used to find an analytical expression for this day. 

T able ( 6 ): The Solar Radiation for 21 st of December [7], 

Time (hour) 



0123456789111111111122222 

012345678901234 

7 

Solar Radiation Intensity ( W/m”) 

0 0 0 0 0 0 0 1 5 1 2 T I 1 T~ 7 1 I T 0 0 0 0 0 0 

837459020288 
5 5 8 3 5 2 2 



It is clear that it is beclouded day, and the sun is shining intermittently, so by using the Fourier series 
an analytical expression can be found and used in any engineering application. 

Figure (5) shows the solar radiation intensity for (27 lh of February) in its analytical expression and 
Meteorological database for (n = 8). 




Figure ( 5): the Meteorological data for solar radiation intensity for (27 th of Feb.) and the analytical 

expression for ( n = 8 ). 

This figure shows the coincidence between them for time period of ( t = 6 to t = 18 ) isn't too close, and 
the coincidence doesn't gives a precise fitting between them. 

Figure (6) shows the solar radiation intensity for the same day but for other values of (n) for (n = 6 and 
n = 10), it is clear from the figure that the plot for (n = 10) is more precision than the plot for (n = 6). 

As a conclusion it is obvious that as ( n ) increases the coincidence precision increases to obtain a good fitness 
between the values for both the tabulated values and the analytical expression. 
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♦ Meteorological 
data 

Fourier (a = 8) 



‘Fourier (a = 

10 ) 



Figure ( 6 ): the Meteorological data for solar radiation intensity for (27 th of Feb.) and the analytical 

expression for ( n = 6 & n = 10 ). 

The plots showing the accuracy with which the (n terms) of this series increase for the given function 
and they are almost indistinguishable from the graph of f(t). 
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